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1.  Introduction 

Many  space  experiments  and  applications  involve 
plasma  beam  emissions.  Electron  beam  experiments 
and  plasma  contactor  experiments  have  been  used 
in  several  ionospheric  active  experiments  to  study 
beam-plasma  interactions,  neutralization  processes,  and 
spacecraft  charging[l,  14,  9],  Ion  beam  experiments 
have  been  conducted  on  SCATHA  to  study  spacecraft 
charging  and  discharging  in  geosynchronous  orbit [3,  4]. 
Electric  propulsion  devices,  such  as  ion  thrusters,  are 
currently  planned  to  use  for  both  commercial  satellites 
and  interplanetary  space  explorations. 

The  purpose  of  this  paper  is  to  study  spacecraft  inter¬ 
actions  induced  by  positive  ion  beam  emissions.  Most 
ion  beam  devices  to-date  are  similar  to  ion  thrusters. 
An  ion  thruster  is  shown  schematically  in  Figure  1.  A 
neutral  gas,  typically  xenon,  is  injected  into  an  ioniza¬ 
tion  chamber,  in  which  the  gas  is  partially  ionized.  The 
ions  are  accelerated  electrostatically  by  a  biased  grid  to 
form  a  high  velocity  beam,  typically  in  the  keV  energy 
range.  Electrons  are  emitted  from  a  neutralizer  to  neu¬ 
tralize  the  ion  beam.  Since  the  ionization  percentage 
is  typically  about  1  to  6%  only,  a  significant  amount  of 
unionized  neutral  gas  will  also  escape  through  the  ion 
optics  and  from  the  neutralizer.  Charge-exchange  col¬ 
lisions  will  occur  between  the  fast  beam  ions  and  slow 
neutrals. 

Ion  beam-spacecraft  interactions  may  be  character¬ 
ized  by  the  following  parameters,  nw/n oo,  £&o/Te, 
Ew/\e$w\,  and  dah/R9C ,  where  n&o  and  £*60  are  the 
beam  ion  density  and  energy  (in  unit  of  eV)  at  beam 
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exit  respectively,  is  the  ambient  plasma  density,  Te 
the  ambient  plasma  temperature  (in  unit  of  eV), 
the  spacecraft  potential,  dah  the  characteristic  sheath 
thickness  surrounding  the  spacecraft,  and  Rac  the  char¬ 
acteristic  dimension  of  the  spacecraft. 

Previous  studies  on  ion  beam-spacecraft  interactions 
were  mostly  performed  for  ion  thrusters  [5,  6,  10, 11, 12, 
18,  20].  During  ion  thruster  operations,  spacecraft  po¬ 
tential  is  typically  kept  low  by  electron  emission  and 
the  beam  ion  energy  is  much  larger  than  the  space¬ 
craft  potential,  E^ o/e$w  >>  1.  Hence,  the  primary 
beam  ions  do  not  return  to  the  spacecraft.  The  ma¬ 
jor  interaction  issue  comes  from  the  charge-exchange 
ions.  The  charge-exchange  ions  can  leave  the  primary 
plume  and  backflow  toward  the  spacecraft.  In  addition 
to  contaminations,  the  charge-exchange  ions  also  rep¬ 
resents  additional  charging  mechanisms,  which  may  in¬ 
duce  plasma  interactions  with  solar  arrays.  Ion  thruster 
experiments  have  been  performed  mostly  either  in  a  lab¬ 
oratory  experiments  or  in  the  low  earth  orbit  environ¬ 
ment,  where  the  ambient  plasma  density  is  relatively 
high.  Other  than  the  thruster  exit  area,  the  spacecraft 
is  typically  covered  only  by  a  thin  sheath  due  to  its 
low  potential,  dah/Rac  «  L  Hence,  our  current  un¬ 
derstanding  of  ion  beam-spacecraft  interactions  mostly 
come  from  an  interaction  regime  that  can  be  character¬ 
ized  by  low-charging  and  thin-sheath  (EbQ/\e$w\  «  1, 
dsh/Rsc  <  1)*  Under  the  low-charging  and  thin-sheath 
condition,  although  the  ion  beam  induced  spacecraft  in¬ 
teraction  is  complex  in  nature,  its  effects  are  expected 
to  be  minimum  from  the  spacecraft  charging  point  of 
view. 

In  this  paper,  we  study  a  different  interaction  regime. 
We  consider  the  situation  that  the  spacecraft  is  very 
negatively  charged,  the  charging  potential  is  compa- 
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rable  or  greater  than  the  initial  ion  beam  energy 
(£W|e*«|  ~  1  or  Ebo/\e$w\  >  1),  and  the  ambient 
plasma  density  is  low  so  that  the  spacecraft  is  coverd 
by  a  thick  sheath,  dsh/R»c  >  1*  We  shall  call  this  in¬ 
teraction  regime  the  high-charging  thick-sheath  regime 
{Ebo/ \e$w\  >  1,  d9h/Rtc  >  1).  Under  the  high-charging 
thick-sheath  condition,  significant  interactions  will  oc¬ 
cur  between  spacecraft  and  both  the  beam  and  charge- 
exchange  ions  within  spacecraft’s  sheath.  Hence,  the  in¬ 
teraction  will  greatly  affect  the  charging  and  discharge 
of  a  spacecraft. 

The  high-charging  condition  can  occur  to  spacecraft 
near  geosynchronous  altitudes,  where  spacecraft  can  be 
charged  negatively  to  the  kV  range  in  the  night  side  or 
during  geomagnetic  storms[7,  8],  A  very  negative  charg¬ 
ing  potential  can  also  occur  when  the  ion  beam  source  is 
operated  with  the  neutralizer  off.  An  example  for  this  is 
the  ion  beam  experiment  on  the  SCATHA  satellite[13]. 
The  thick-sheath  condition  occurs  to  spacecraft  in  the 
geosynchronous  orbit  plasma  environment  or  interplan¬ 
etary  plasma  environment,  where  a  spacecraft  is  always 
covered  by  a  thick  sheath  due  to  the  low  ambient  den¬ 
sity. 

This  study  is  based  on  fully  3-dimensional  (3-D)  com¬ 
puter  particle  simulations.  Recently,  Samanta  Roy  et  al. 
used  electrostatic  Particle-in-Cell  (PIC)  simulations  to 
model  the  far- downstream  region  of  an  ion  thruster  and 
study  charge  exchange  ion  backflow[17,  18].  Wang  et 
al.  developed  full  particle  Particle-in-Cell  with  Monte 
Carlo  collision  (PIC-MCC)  models  for  single  thruster 
as  well  as  multiple  thruster  plasma  plumes  under  both 
ground  test  and  in-space  conditions [20,  21],  The  em¬ 
phasis  of  these  studies  was  to  characterize  the  plasma 
plume  and  the  charge  exchange  ions  in  the  down  stream 
region.  Here,  we  move  towards  a  more  detailed  model  of 
spacecraft  interactions  induced  by  ion  beam  emission  by 
considering  ion  beam  emissions  from  a  3-dimensional,  fi¬ 
nite  size,  negatively  charged  spacecraft  emersed  in  space 
plasmas.  The  emphasis  is  to  study  high-charging  thick- 
sheath  interactions  induced  by  ion  emission  without  neu¬ 
tralizing  electrons.  We  discuss  the  plasma  environment 
surrounding  the  entire  spacecraft  and  the  backflow  of 
beam  and  charging- exchang  ions  under  various  space¬ 
craft  charging  conditions. 

2.  Formulation  and  Approach 

Our  simulation  setup  is  shown  in  Fig.2.  The  space¬ 
craft  is  modeled  as  a  3-dimensional  box  structure  with 
a  conducting  surface  and  a  surface  potential  relative 
to  the  ambient.  We  assume  that  the  ambient  plasma  is 
a  H+  plasma  represented  by  a  single  Maxwellian  distri¬ 
bution.  (rioo,  Ti,  and  Te  denote  the  ambient  plasma 


density,  ion*  and  electron  temperature  respectively.) 
The  ion  beam  source  is  taken  to  be  a  30cm  Xenon  ion 
thruster,  (n&o  and  t;*)  denote  the  average  beam  ion  den¬ 
sity  and  beam  velocity  at  thruster  exit  respectively.) 

For  convenience  of  the  analysis,  we  normalize  the 
spatial  dimension  by  the  Debye  length  at  thruster  exit 
Aoo  =  yjKTJfambQe2,  time  by  the  xenon  ion  plasma 
frequency  at  thruster  exit  upiQ  =  y/iirnto  /mxenon ,  ion 
velocity  by  the  xenon  ion  acoustic  velocity  within  the 
ion  beam  C9  =  y/Te/mxenon,  and  the  potential  by  the 
ambient  electron  temperature  Te. 

The  equilibrium  charging  state  for  a  uniformly  con¬ 
ducting  body  is  determined  by  the  total  current  balance 
condition:  HI(<&W)  —  0.  For  our  problem,  the  charging 
equation  is  written  as 

^f(^tu)  ~  I\  —  Ie  ~~  Itri  “b  f&e  ~b  Ireturn  ~  0  (1) 

where  and  Ic  are  ion  and  electron  current  from  the 
ambient  plasma  respectively,  is  the  beam  ion  current, 
and  7^  is  the  electron  current  emitted  from  the  neutral¬ 
izer.  Iretum  is  the  beam  induced  return  current,  which 
includes  both  the  beam  ion  current  hw  and  the  charge 
exchange  ion  current  Icexr  collected  by  the  spacecraft. 

From  the  probe  theory,  the  ambient  electron  and  ion 
current  to  a  negatively  charged,  stationary  conducting 
sphere  under  the  thick-sheath  condition  is  given  by[15]: 

Ic  —  J eoo  exp{$w/Te)A,e  =  en°o\jj~~  exp($w/Te)A,c 

(2) 

It  =  JiooO-  -  $w/Ti)A,c  -  en^  -  $w/Ti)A,c 

(3) 

where  A8C  is  the  spacecraft  area.  However,  to  calculate 
Ireturn,  one  needs  to  obtain  the  self-consistent  solution 
of  space  charge  distribution,  electric  field,  and  particle 
orbits. 

We  have  developed  a  set  of  full  particle  and  hybrid 
PIC-MCC  simulation  codes  to  simulate  ion  thruster 
plasma  interactions  [20].  The  numerical  models  and  al¬ 
gorithms  have  been  discussed  in  detail  in  Ref[20,  22,  23]. 
Since  the  focus  of  this  paper  is  on  ion  interactions  with 
a  negatively  charged  body,  simulation  results  presented 
in  this  paper  are  derived  from  a  hybrid  PIC-MCC  code. 
The  electron  density  is  assumed  to  follow  the  Boltz¬ 
mann  distribution  nc  =  nrcf  exp(($  -  $rc/)/Tc),  where 
nre/  and  $re/  are  the  electron  density  and  potential  at 
a  reference  point  respectively.  In  particular,  the  ambi¬ 
ent  electron  density  surrounding  a  negatively  charged 
spacecraft  is  given  by 

ne  =  exp  (e$/KTc)  (4) 
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The  ambient  ion  density,  which  is  much  lower  than  the 
beam  ion  density,  is  assumed  to  be  a  constant  back¬ 
ground.  The  beam  ions  and  charge-exchange  ions  are 
treated  as  test  particles. 

At  each  time  step,  the  beam  ions  are  injected  into  the 
simulation  domain  from  the  beam  source  exit.  Experi¬ 
ments  have  shown  that  the  beam  ions  are  emitted  with 
a  spatial  density  distribution  close  to  that  of  a  Gaus¬ 
sian  distribution  and  a  divergence  angle  of  about  15°  to 
20°.  To  simulate  the  ion  beam  current  emitted  from  a 
thruster,  the  particles  are  injected  in  such  a  way  that 
the  resulting  flux  has  a  density  in  Gaussian  distribution 
at  the  exit  plane 

|Ai|  =  JbmaXexp(-(r/rT)2),  r  <  rT  (5) 

where  Jbmax  is  the  density  at  the  thruster  exit  center,  r 
is  the  distance  to  the  center  on  the  thruster  exit  plane, 
and  rr  is  the  thruster  radius.  In  this  study,  we  take  the 
divergence  angle  to  be  15°. 

In  an  ion  thruster,  not  all  the  propellant  is  ionized  due 
to  discharge  chamber  inefficiencies  and  operating  con¬ 
straints.  Along  with  the  beam  ions,  significant  amount 
of  neutrals  also  flow  out  of  the  thruster  and  form  a  neu¬ 
tral  plume.  In  the  simulation,  the  neutral  plume  is  mod¬ 
eled  as  a  steady  background.  As  in  [17],  the  density 
distribution  of  the  neutral  plume  is  modeled  as  that  of 
a  free  molecular  flow  from  a  point  source  located  at  one 
thruster  radius  tt  behind  the  thruster  exit [17].  Hence, 
the  density  distribution  of  the  neutral  plume  is  given 
by: 

nn(R, 6)  =  ann0(l  -  (1  +  (^)2)-1/2)  cos  6  (6) 

XL 

where  nn o  is  the  average  neutral  density  at  thruster  exit, 
R  is  the  distance  to  the  point  source,  9  is  the  angle  be¬ 
tween  R  and  the  downstream  axis,  and  a  is  a  correction 
factor.  UnQ  is  related  to  the  ion  beam  via  the  following 
relationship  [17]: 

47160^60,1  -v 

n„0  -  — tt- . (~  "  1)  (7) 

C  Tip 

where  C  =  y/ &kTw  / 7rmxenon  is  the  mean  thermal  speed 
of  the  neutrals,  and  t]p  is  the  propellant  utilization  effi¬ 
ciency.  Typically,  the  thruster  discharge  chamber  wall 
temperature  is  Tw  ~  500K.  For  normal  thruster  oper¬ 
ations,  Tjp  range  from  0.7  to  0.9.  A  significant  higher 
neutral  plume  density  may  also  be  obtained  either  due 
to  low  propellant  utilization  or  through  extra  neutral 
sources. 

Charge-exchange  collisions  will  occur  between  beam 
ions  and  neutrals.  The  charge-exchange  collision  cross 


section  for  xenon  can  be  expressed  by[16,  17]: 

acex  =  (*h  In  i>w  +  fc*)2  x  10”20  m2  (8) 

where  vu  is  beam  ion  velocity  with  unit  in  m/s,  k\  = 
-0.8821,  and  k2  =  15.1262. 

As  in  a  standard  PIC-MCC  code[2],  the  trajectory  of 
each  particle  is  integrated  from 

f  f  =  f  (9) 

Citr  C  dir 

using  a  standard  leapfrog  scheme.  The  probability  that 
a  charged  particle  suffers  a  collision  within  time  t  is 
given  by 

P(t)  —  l  —  exp(-  f  v(t)dt)  (10) 

Jo 

where  v  —  nn{x)va(y)  is  the  collision  frequency,  and 
a(v)  is  given  by  eq(8).  Since  the  neutral  density,  which 
is  defined  on  grid  points,  is  nonuniform,  the  collision 
frequency  for  each  particle  is  obtained  by  interpolating 
the  neutral  density  nn(z,y,z)  to  the  particle  position, 
similar  to  the  field  interpolation  in  a  PIC  code.  At  each 
time  step,  for  each  particle,  the  accumulated  collision 
probability  in  the  time  step  is  calculated,  and  a  random 
number  PTan  evenly  distributed  between  0  and  1  is  then 
chosen  to  determine  whether  a  collision  has  occured. 
The  self-consistent  electric  field  is  than  obtained  from 
the  Poisson’s  equation 

V2$  =  — 47rp  (11) 

Finally,  since  we  are  simulating  a  thick  sheath  situa¬ 
tion,  computationally  it  is  not  feasible  to  set  the  simu¬ 
lation  domain  large  enough  for  the  outer  boundary  to 
be  the  undisturbed  ambient.  Hence,  a  Neumann  con¬ 
dition  V^n  =  0  is  used  at  all  outer  boundaries  of  the 
simulation  domain.  In  this  paper,  the  number  of  grid 
points  used  for  the  simulation  domain  is  47  x  47  x  47 
and  the  number  of  test  particles  is  in  the  range  of  106. 

3.  Results  and  Discussions 

In  this  section,  we  discuss  simulation  results.  We  con¬ 
sider  the  spacecraft  configuration  shown  in  Figure  2a 
(no  solar  array).  The  spacecraft  is  taken  to  be  a  cubic 
box  with  dimension  lm  x  lm  x  lm.  The  ion  beam  source 
is  taken  to  be  a  30cm  Xenon  ion  thruster.  The  beam 
ions  exit  the  thruster  with  an  energy  of  E^o  =  800eV 
Since  our  main  purpose  here  is  to  study  ion  beam  in¬ 
teractions,  the  ion  beam  current  is  taken  to  be  much 
less  than  the  value  used  for  propulsion  applications. 
To  cover  a  range  of  interaction  parameters,  ri6o/noo) 
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EM/Tei  £W|e$u;|,  and  dah/Rac,  we  present  four  simu¬ 
lation  cases,  as  summarizd  in  Table  1. 

In  simulation  cases  A  and  B,  we  take  the  ambient 
plasma  environment  to  be  a  “cold”  plasma  with  ~ 
102  cm"3,  Ti  =  0.5 Tc  and  Tc  ~  8eV  ~  100). 

(Ambient  current  densities:  ~  1.3xl0”7  A/m ?  and 

Jeoo  —  7.5  x  10”6  A/m2.)  We  take  the  average  beam  ion 
density  at  thruster  exit  to  be  n&o  ~  103noo  ~  105  cm”3. 
(Average  beam  current  density  at  exit  Jbo  =  |e|n6ov&o  ~ 
5.3  x  10“4  A/m2.)  In  case  A,  we  assume  the  ion  beam 
source  is  operated  with  the  neutralizer  on,  and  hence, 
the  beam  is  quasineutral.  Since  there  is  no  charging 
mechanism  involved  in  this  situation,  the  spacecraft  can 
be  charged  only  to  a  low  potential.  We  take  the  space¬ 
craft  potential  to  be  -$W/TC  ~  3  (2?M)/|e$w|  ~  33.3). 
Hence,  case  A  is  a  low-charging  case. 

In  case  B,  we  assume  that  the  ion  beam  source  is  op¬ 
erated  with  the  neutralizer  off.  Hence,  the  ion  beam  is 
only  neutralized  by  the  ambient  electrons.  (The  ambi¬ 
ent  electron  density  is  given  by  eq(4).)  It’s  interesting 
to  estimate  the  maximum  charging  potential  a  space¬ 
craft  can  reach  under  this  condition.  From  eq(l),  the 
maximum  charging  potential  is  given  by  —  hi  —  0,  as¬ 
suming  there  are  no  returning  ion  current  and  electron 
current.  Hence,  we  have 

<  1  7ifro  At  I  mp  FTj  [2Ebo  ^ 

Te  27T  tIqq  Aac  y  mxenon  y  Te\/  Te 

where  At  is  the  ion  beam  emitting  area.  For  the  pa¬ 
rameters  used  here,  we  found  $W/TC - 47.  We  shall 

use  this  value  for  in  case  B  {Ebo/\e$w\  ~  2.1). 

In  simulation  cases  C  and  D,  we  take  the  ambient 
plasma  environment  to  be  a  “warm”  plasma  with  ~ 
102  cm"3,  Ti  =  0.5TC  and  Te  ~  80eV  (£W/TC  ~  10). 
(Ambient  current  densities:  Jioo  ~~  4,1  x  10"7  A/m 2 
and  Jcoo  ~  2.4  x  10”5  A/m2.)  We  take  the  aver¬ 
age  beam  ion  density  at  thruster  exit  to  be  — 
104noo  ~  106  cm”3.  (Average  beam  current  density 
at  exit  J&o  —  5.3  x  10”3  A/m2.)  Again,  we  assume  that 
the  ion  beam  source  is  operated  with  the  neutralizer 
off,  and  the  spacecraft  is  charged  to  some  very  negative 
potential.  In  the  high-charging  case,  it  is  not  possible 
to  estimate  the  charging  potential  a  prio  because  the 
return  current  can  no  longer  be  neglected.  In  case  C, 
we  take  —$W/Te  ~  10,  i.e.  Ebo/\e$w\  ~  1.  In  case  D, 
we  take  ~$W/Te  ~  75,  i.e.  Eb0/\e$w\  ~  0.13. 

Since  the  ambient  Debye  length  is  larger  than  the 
spacecraft  dimension,  the  spacecraft  is  covered  by  a 
thick  sheath  in  all  simulation  cases  here.  In  the  simula¬ 
tion  results  presented  here  (with  the  exception  of  Figure 
3),  the  spacecraft  is  located  at  2  <  x  <  16,  15  <  y  <  29, 


Cases 

Ebo/\e$w\ 

Tlbo/Hoo 

Em/T' 

-*w/Te 

A 

33.3 

103 

100 

3 

B 

2.1 

103 

100 

47 

C 

1.0 

104 

10 

10 

D 

0.13 

104 

10 

75 

Table  1:  Simulation  cases 


and  15  <  z  <  29.  The  thruster  exit  center  is  located  at 
x  =  18,  y  =  22,  and  z  =  22.  The  grid  resolution  is  given 
by  the  Debye  length  at  thruster  exit  \do  ~  6.7cm. 

The  simulation  results  are  shown  from  Figure  3 
through  Figure  12.  In  the  following,  most  contour  plots 
are  for  a  xy  plane  cutting  through  the  spacecraft  and 
thruster  center  at  z  =  Zthmster  =  22.  Hereafter,  this 
plane  will  be  referred  to  as  the  “center  xy  plane”.  All 
particle  plots  are  for  particles  located  within  a  layer  of 
±1  cell  of  the  center  xy  plane.  Hereafter,  this  layer  will 
be  referred  to  as  the  “center  layer” . 

Figure  3  shows  the  potential  contours  on  the  center 
xy  plane  for  the  no  beam  situation.  Figure  4  shows  the 
neutral  plume  density  contours  on  the  center  xy  plane. 
Since  the  neutral  density  decreases  as  J?”2,  most  of  the 
charge-exchange  ions  are  generated  close  to  the  thruster 
exit. 

The  results  from  simulation  case  A  are  shown  in  Fig¬ 
ure  5  ($  contours  and  £  field  vectors  on  the  center  xy 
plane)  and  Figure  6  (beam  ions  and  charge-exchange 
ions  within  the  center  layer).  This  case  is  similar  to 
the  ion  thruster  situation.  Since  Ebo/\e$w\  »  1  and 
the  beam  is  quasineutral,  the  beam  ions  are  not  influ¬ 
enced  by  the  potential  field.  The  structure  of  the  poten¬ 
tial  field  is  dominated  by  the  disturbance  from  the  high 
density  ion  beam.  The  potential  in  the  beam  center 
$  center  and  the  ambient  potential  has  the  following 
approximate  relation: 

$  center  ~  ^oo  —  relnnjeom/n  oo 

Further  away  from  the  thruster,  the  potential  diffuse  to 
that  of  the  ambient.  Collisions  between  the  beam  ions 
and  the  neutral  background  generate  charge  exchange 
ions.  The  charge  exchange  ion  production  rate  is  pro¬ 
portional  to  the  neutral  density  shown  in  Figure  4.  In 
contrast  to  the  beam  ions,  the  motion  of  the  charge  ex¬ 
change  ions  are  influenced  by  the  potential  field  due  to 
their  low  kinetic  energy.  The  potential  field  influences 
the  charge  exchange  ions  in  two  ways.  First,  charge 
exchange  ions  produced  near  the  thruster  exit  will  be 
accelerated  towards  the  thruster  exit.  This  backflow 
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may  cause  potential  erosion  on  the  acceleration  grids. 
Second,  since  the  plume  center  has  a  higher  potential, 
charge  exchange  ions  produced  within  the  beam  can  flow 
radially  outward  the  beam  region.  It  is  well  recognized 
that,  once  outside  the  plume,  charge  exchange  ions  may 
become  a  potential  contamination  source.  Except  for 
those  charge- exchange  ions  generated  right  downstream 
of  the  thruster,  the  collection  of  charge-exchange  ions 
by  the  spacecraft  is  similar  to  the  current  collection  by 
a  probe  in  the  orbital  motion  limited  regime.  As  the  E 
vector  plot  in  Figure  5b  and  charge-exchange  ion  plot 
in  Figure  6b  clearly  indicate,  a  symmetrically  shaped 
spacecraft  charged  to  a  few  Te  can  not  collect  a  signifi¬ 
cant  amount  of  charge-exchange  ion  current. 

The  results  from  simulation  case  B  are  shown  in  Fig¬ 
ure  7  (<£  contours  and  E  field  vectors)  and  Figure  8 
(beam  ions  and  charge- exchange  ions  within  the  center 
layer).  Figure  7  clearly  shows  that  the  structure  of  the 
potential  field  is  a  combined  result  of  the  surface  poten¬ 
tial  and  the  beam  space  charge.  The  surface  potential 
generates  a  sheath  covering  the  spacecraft.  On  the  other 
hand,  the  potential  structure  generated  by  the  beam  is 
similar  to  the  expansion  fan  structure  associated  with 
a  plasma  expanding  into  a  vacuum[19].  The  transition 
of  the  potential  structure  from  that  dominated  by  the 
sheath  to  that  dominated  by  the  beam  is  evident  in 
Fig.  7.  Since  le^l  and  o  is  comparable  in  this  case 
(EbQ/\e$w\  ~  2.1),  the  beam  ions  are  affect  by  the  po¬ 
tential  field  (the  slowing  down  of  beam  ions  are  evident 
in  Figure  8c).  The  E  field  structure  shown  in  Figure  7b 
and  charge  exchange  ion  plot  in  Figure  8  also  indicate 
that  the  charge-exchange  ions  not  only  flow  radially  out¬ 
ward  the  beam  region  but  also  backflow  to  orbit  around 
the  spacecraft.  The  charg-exchange  ions  collected  by 
the  spacecraft  is  ICexr/Iioo  ~  3  x  10~7. 

The  results  from  simulation  case  C  are  shown  in  Fig¬ 
ure  9  and  Figure  10.  Due  to  the  high  beam  density,  the 
sheath  near  the  thruster  exit  is  much  thinner  than  the 
sheath  covering  the  spacecraft.  Since  Et,o/\e$w\  ~  1 
in  this  case,  the  beam  ions  are  stoped  and  pulled  back 
over  the  distance  of  the  thruster  sheath.  Hence,  in  ad¬ 
dition  to  the  charge-exchange  ions,  the  spacecraft  also 
collects  the  beam  ions.  In  case  D,  we  further  increase 
the  spacecraft  potential  so  that  Ebo/\e$w\  ~  0.13  (Fig¬ 
ures  11  and  12).  Case  D  is  the  opposite  limit  of  case 
A  in  that  the  structure  of  the  potential  field  is  dom¬ 
inated  by  the  spacecraft  sheath.  The  ion  beam  only 
represents  a  small  perturbation.  The  entire  population 
of  the  beam  ions  and  charge-exchange  ions  sure  trapped 
within  the  potential  well  of  the  spacecraft  and  will  be 
collected.  Of  course,  this  will  led  to  the  discharging  of 
spacecraft. 


4.  Summary  and  Conclusions 

We  have  developed  a  3-D  PIC-MCC  simulation  model 
to  study  ion  beam-spacecraft  interactions.  Simulation 
results  are  presented  for  ion  beam  emissions  under  sev¬ 
eral  charging  conditions.  The  interaction  characteristics 
are  discussed  for  both  the  high-charging  and  the  low- 
charging  regime.  When  the  interaction  is  within  the 
low-charing  regime  (Ebo/\e$w\  »  1),  such  as  that  in 
an  ion  propulsion  application,  the  interactions  mainly 
come  from  the  charge-exchange  ions.  We  find  that 
charge-exchange  ions  typically  does  not  alter  a  space¬ 
craft’s  charging  environment.  When  the  interaction  is 
within  the  high-charing  regime  {Ebo/\e$w\  <  1),  both 
the  charg-exchange  ions  and  the  beam  ions  will  return 
to  the  spacecraft.  The  returning  beam  ions  can  led  to 
the  discharging  of  the  spacecraft.  It  is  also  interesting 
to  note  that  ion  beam  emissions  in  the  high-charging 
regime  may  be  used  to  discharge  a  spacecraft  efficiently. 
For  instance,  a  spacecraft  often  have  “hot”  spots  on 
dielectric  surfaces  and  it  is  difficult  to  discharge  these 
spots  by  electron  emission  from  the  metallic  ground.  In 
this  situation,  not  only  the  ions  return  to  the  spacecraft 
but  they  tend  to  land  on  these  “hot”  spots,  where  the 
potential  is  more  negative.  This  reduces  the  level  of 
differential  charging  gradually.  While  this  study  has  fo¬ 
cused  on  the  returning  current  under  a  given  spacecraft 
potential,  in  future  we  plan  to  include  a  dynamic  cal¬ 
culation  of  the  spacecraft  potential  and  study  the  time 
development  of  spacecraft  charging/discharge  during  an 
ion  beam  emission.  Other  interaction  issues,  such  as 
plasma  waves  and  instabilities  induced  by  an  ion  beam, 
will  also  be  studied. 
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Figure  Captions 


Figure  1:  Ion  thruster  (Upper).  Space  charge  poten¬ 
tial  along  the  ion  beam  (Lower). 

Figure  2:  Model  setup,  a)  spacecraft,  b)  spacecraft 
with  solar  arrays. 

Figure  3:  Potential  contours  on  the  center  xy  plane 
for  no  beam  situation. 

Figure  4:  Neutral  plume  density  contours  on  the  cen¬ 
ter  xy  plane. 

Figure  5:  Potential  contours  (a)  and  E  field  (b)  on 
the  center  xy  plane.  Simulation  case  A. 

Figure  6:  Beam  ion  positions  (a)  and  charge-exchange 
ion  positions  (b)  within  the  center  layer.  Simulation 
case  A. 

Figure  7:  Potential  contours  (a)  and  3  field  (b)  on 
the  center  xy  plane;  Potential  contours  on  the  z  =  15 
xy  cutting  plane  (c).  Simulation  case  B. 

Figure  8:  Beam  ion  positions  (a),  charge-exchange 
ion  positions  (b),  and  beam/charge-exchange  ion  vx  vs. 
x  (c)  within  the  center  layer.  Simulation  case  B. 

Figure  9:  Potential  contours  on  the  center  xy  plane 
(a)  and  on  the  z  =  15  xy  cutting  plane  (b).  Simulation 
case  C. 

Figure  10:  Beam  ion  positions  (a),  beam  ion  vx  vs. 
x  (b),  charge-exchange  ion  positions  (c),  and  charge- 
exchange  ion  vx  vs.  x  (d)  within  the  center  layer.  Sim¬ 
ulation  case  C. 

Figure  11:  Potential  contours  on  the  center  xy  plane 
(a)  and  on  the  z  =  15  xy  cutting  plane  (b).  Simulation 
case  D. 

Figure  12:  Beam  ion  positions  (a),  beam  ion  vx  vs. 
x  (b),  charge-exchange  ion  positions  (c),  and  charge- 
exchange  ion  vx  vs.  x  (d)  within  the  center  layer.  Sim¬ 
ulation  case  D. 
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a)  $  Cutting  plane  at  ZVIEW=  22 


a)  beam  ion  position 


b)  E  Cutting  plane  at  ZVIEW=22. 


b)  charge-exchange  ion  position 


20.  x  40. 


c)  $  Cutting  plane  at  ZVIEW=  15. 


c)  phase  space 


beam  ions 


charge- exmange  ions 


,)  4>  Cutting  plane  at  ZVIEW=  22. 


b)  <J>  Cutting  plane  at  ZVIEW=  15 


